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We present a new hybrid method to solve the relativistic Hartree-Fock-Roothan equations where 
the one- and two-electron radial integrals are evaluated numerically by defining the basis functions 
00 ■ on a grid. This procedure reduces the computational costs in the evaluation of two-electron radial 

ON ' integrals. The orbitals generated by this method are employed to compute the ionization potentials, 

excitation energies and oscillator strengths of alkali-metal atoms and elements of group IIIA through 
second order many-body perturbation theory (MBPT). The computed properties are in excellent 
O ^ agreement with the experiment and other correlated theories. 
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^ ■ I. INTRODUCTION 

The two most critical choices in the application of many-body perturbation theory (MBPT) (relativistic as 

well as non-relativistic) to atomic and molecular systems involve the appropriate selection of basis functions and 
the partitioning of the full Hamiltonian H into a zeroth order Hamiltonian H and a perturbation V These 
choices become extremely important when highly accurate estimates of various properties are demanded from low 
<Z2 ■ order perturbative computations. Intensive research has focused on developing appropriate basis sets and methods 
to minimize the error between the theoretically computed properties and its experimental value. The strong 
dependence of the convergence of MBPT on the choice of Hq was first demonstrated by Kelly [Q in his pioneering work 
^**' on beryllium atom. Using a V N ~ X instead of the traditional V N potential for the excited orbitals, Kelly obtained a 
vast improvement in the perturbative convergence for that atom. He also demonstrated that more rapid convergence 
can be achieved from a shifted denominator that corresponds to the summation of a certain class of diagrams to all 
order. 

It is well-known that the theoretical treatment of the heavy atoms must incorporate certain special features that 
\ are not essential for the light atoms. This is largely due to the fact that the relativistic effects are so large for heavy 
atoms that it is imperative to treat them by using the relativistic Dirac equation. Despite its enormous computational 
complexity and cost, tremendous progress has been made over the past decade and a half in solving the four-component 
Q\ ■ Dirac equations for many-electron systems using numerical Dirac-Fock (DF) and the finite basis set expansion (FBSE) 
method (9 - 21 . The numerical atomic DF self-consistent field (SCF) calculation is more compact and accurate but 
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its extension to molecular systems (multi-center many-electron systems) is cumbersome. Moreover, the generation of 
virtual orbital is tedious and frequently encounters convergence difficulties. The FBSE method, on the other hand, is 
C/3 | rather simple and its extension to molecules is straightforward. Also, the generation of occupied and virtual orbitals 
. do not require separate computations. 

The success of the relativistic FBSE method lies in its proper imposition of kinetic-balance condition |^2j between 
the large and small component spinor which in essence can be regarded as a proper boundary condition upon the basis 
set. Several papers by Grant et al. [ fi"3[ and Parpia et al |23| ] among others document the success of the relativistic 
FBSE method. However, in their finite basis set calculation for light to heavy atoms (Z=2-50 and 80), Grant et 
al. jl3| employed kinetically balanced Slater-type orbitals (STOs) which have the correct functional behavior but 
are particularly unsuitable for analytical self-consistent field (SCF) molecular calculations. Gaussian type orbitals 
(GTOs) or contracted Gaussian type orbitals (CGTO), on the other hand, are suitable in the evaluation of multi- 
center integrals in molecules. It was shown by Ishikawa et al. that GTOs can give rise to a natural description of 
the relativistic wave- functions within a finite nucleus. 

The most important feature of the FBSE method in STO [exp(— £r)] or GTO [exp(— (r 2 )] framework is to determine 
the appropriate exponential parameter because the quality of the wave- function largely depends upon this parameter 
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and in recent years there has been an increased interest in finding out the appropriate exponential parameter and 
contraction coefficients (for CGTO) that can provide correct functional behavior of the relativistic wave- functions at 
the nucleus Jl5| , p5| -[2^] . For instance, Matsuoka et al. |^6| have reported accurate configuration average DF energies 
for various atoms through the FBSE method using kinetically balanced well-tempered basis set in the framework 
finite nuclear size approximation. While Matsuoka et al. p6| ] used a well-tempered Gaussian basis set in computing 
DF energies, Clementi et al. flEf employed kinetically balanced geometric-type exponent for the Gaussian primitives 
and obtained DF energies for various atoms that are comparable to the numerical DF value [^9). Later Malli et 
al. |^] reported all-electron ab-initio fully relativistic DF and DF-Breit calculations for polyatomic systems using a 
relativistic universal Gaussian basis set and recently Pernpointer et al. |27| employed a relativistic CGTO basis set 
in their relativistic coupled cluster calculation for the nuclear quadrupole moment of CsF. Though the slater (STO) 
and Gaussian (GTO) types of basis functions are most widely used in atomic many-body calculations, this choice is, 
in principle, arbitrary. Since, it is beyond the scope of this present work to discuss this aspect at length, we refer the 
review articles by Grant |30j and Sapirstein |3l|] for details. 

We have developed a numerical procedure to solve the atomic relativistic DF-SCF equations using the FBSE method. 
This new approach is basically a hybrid of numerical and analytical DF (FBSE) methods. Here, the large and small 
component radial functions are expanded in terms of Gaussian primitives on a grid using appropriate constraints 
on the small component radial basis to impose the kinetic-balance condition. While the large and small component 
part of the radial functions are generated (on a grid) through the FBSE procedure, the one- and two-electron radial 
integrals are evaluated numerically to avoid the complicated analytical expressions for the two-electron direct and 
exchange radial integrals (the analytical evaluation of one electron radial integral is rather straight forward) . This is 
the part which differs from the conventional FBSE method for solving DF-SCF or HF-SCF equations. This procedure 
(numerical computation of two-electron integrals) also provides an easy way to reduces the N * (N + l)/2 operations 
to N c operations (TV and N c corresponds to the number of basis set and occupied orbitals, respectively) in DF-SCF 
computation, and, thereby reduces the computational time of relativistic self-consistent field calculations for heavy 
atoms. In the perturbative computations of ground and excited state properties, the two-electron radial integrals are 
also directly computed (numerically) wherever they appear to avoid the two-electron integral storage problem. In 
this paper, we present some pilot calculations of the ionization potentials and excitation energies of alkali-metal and 
elements of group IIIA computed through second order MBPT using relativistic wave-functions obtained from the 
hybrid DF-SCF approach. 

In Sec. II, we describe the hybrid DF-SCF method that has been used to generate the relativistic single-particle 
atomic orbitals for post-Dirac-Fock computations. Sec. Ill briefly reviews the background of the MBPT approach for 
computing ionization potentials (IP), electron affinities (EA) and excitation energies (EE). The numerical results are 
presented in Sec. IV and compared with other perturbative calculations where available. We make some concluding 
remarks in Sec. V. 



II. HYBRID RELATIVISTIC HARTREE FOCK ROOTHAN EQUATION 



The Dirac-Coulomb Hamiltonian for a many-electron system can be conveniently written as 



H = V[ca 4 .p, + (ft - l)mc 2 + V nuc {n)} + 1 V £ (2.1) 
in which the Dirac operators a and (3 are expressed by the matrices 

% o) - ,) M 

where a stands for the Pauli matrices and / is the 2x2 unit matrix. 

In the central field approximation, the SCF equations are determined by minimizing the energy functional E with 
respect to $, where E is given by 

N 1 2 

E = ($| Y\ c ^-Pi + (A - l ) m <? + Vnudn)] + - J *) (2-3) 

2—1 l^J ' J ' 

and determinantal wave-function (antisymmetric) u is built from single particle orbitals 
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where r _1 P nK (r) and r _1 Q nK (r) are the large and small component radial wave- functions, respectively that satisfy 
the orthonormality condition 

/>oo 

dr[P nK {r)P n , K {r) + Q nK {r)Q n , K {r)\ = 5 nn , (2.5) 

Here, the quantum number k classifies the orbital according to their symmetry and is given by 

K = -2(j-l)(j+±) (2.6) 

where I is the orbital quantum number and j = i±i is the total angular quantum number. Here, the spinors X« m (0) <j>) 
are given 

X™ = X! ^(4-?' m ~ a ' °) Y Lrn-a{0, <j>)r} a (2.7) 

where C(lhj\m — a, a) and Y^ m _ a (d, cf>) represent the Clebsch- Gordon coefficients and the normalized spherical 
harmonics, respectively, and the r\ a stands for the two-component spinors. 

With these definitions, it can be easily shown that the application of the variation principle to Eq. (2.3) leads to 
a coupled integro-differential equations in P nK (r) and Q nK (r). Therefore, to obtain the numerical wave- functions, we 
have to solve these two coupled integro-differential equations. Alternatively, a pseudo-eigenvalue equation (Hartree- 
Fock-Roothan) |32] | can be obtained by using an analytic expansion-type wave-functions and minimizing the energy 
functionals E with respect to the expansion coefficients. 

It has been found that the numerical wave-functions have more accurate asymptotic behavior than the analytical 
ones, though both provide total energies of comparable accuracy. The accuracy of the total energy and wave-function 
obtained through the Dirac-Fock-Roothan equation (FBSE method) can in principle be enhanced to any degree 
by increasing the number of basis functions, but in reality only a finite number of basis can be used because the 
computational time increases very rapidly with the increasing number of basis functions. Moreover, the use of large 
basis functions severely impedes the efficiency of the post-Dirac-Fock computations. 

In the present paper, we introduce a hybrid scheme to solve the DF equation through the pseudo-eigenvalue approach 
where basis functions are defined on a grid and one- and two-electron radial integrals are evaluated numerically as 
opposed to the conventional relativistic Hartree-Fock-Roothan equations. Since, the basis functions are defined on a 
grid and the matrix elements appearing in the relativistic Hartree-Fock-Roothan equations are evaluated numerically, 
this scheme can be regarded as a combination of numerical and analytical approach to the solution of DF-SCF 
equation. Here, like the traditional analytical basis set expansion approach, the large and small components of the 
radial wave-functions are expressed as linear combination of basis functions, i.e., 

P 

and 

g„ K (r)=^Cf pff f p (r) (2.9) 
p 

where the summation index p runs over the number of basis functions N, g^ p (r) and g^ p (r) are basis functions belonging 
to the large and small components, respectively, and C^ p and C^ p are the corresponding expansion coefficients. 
Though, any basis functions can be used, we have chosen Gaussian-type of orbitals (GTOs) that has the following 
form for the large component, 

ff«p(r)=^r"-e- a "- a (2.10) 

with 

a p = a a f3'P- 1 (2.11) 
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where ao,(3 are user defined constants, n K specifies the orbital symmetry (1 for s, 2 for p, etc.) and A/jf is the 
normalization factor for the large component. The small component part of the basis function is obtained by imposing 
the kinetic balance and has the following form 



»-K S (i + 7^) (2-12) 



' dr 



K S = \Ur JL ^l^ +K + n K )-l] (2.13) 



where 

\ 

" " 2n K - 1 

Using the above definitions, the Dirac-Fock-Roothan equation for closed shell system can be cast into a pseudo- 
eigenvalue equation of the form 

FC = SCe (2.14) 

where F is the Fock matrix and S, C and e are overlap, eigenvector and eigenvalue matrices, respectively. This 
pseudo-eigenvalue equation is first transformed into an eigenvalue equation F'C = C'e, which on diagonalization 
produces the desired eigenvalues (e) and eigenvectors (C = S~ 1 / 2 C). Since, the detailed derivation of relativistic 
Hartree-Fock-Roothan equation and its matrix elements (analytical as well as numerical form) has been discussed in 
details by several authors dJl^l , we conclude this section by reiterating that the DF matrix elements appearing in the 
hybrid relativistic Hartree-Fock- Roothan equations are evaluated numerically to avoid the evaluation of complicated 
analytical expression of two-electron matrix elements and to improve the accuracy of the orbital properties. The 
present procedure also provides an easy route to implement N c (N c specifies the number of occupied orbitals) opera- 
tions instead of iV 2 (N denotes the number of basis functions) for the evaluation of the two-electron radial integrals 
that appear in DF-SCF equation. A brief outline of the scheme is the following: 

In the SCF procedure, the integrals and the matrices are evaluated over the members of the basis set {4>fj,} rather 
than over the members of the set of solutions {4>i} because the atomic or molecular orbitals (solutions of SCF 
equations) are not known until the calculation is complete. Since these two sets of functions are related by 

N 

i>i = J2 C ^ (2.15) 
the two-electron matrix element of F (the Hartree-Fock potential term) in {cf>} basis can be written as 

un = yv&Vci— iMc) =EEE c ^*i-i^> (2.i6) 
c ri2 cir, ri2 

which involves a two-index tranformation. However, this two-index transformation process can be easily avoided 
by evaluating the Uij matrix elements in a mixed basis i.e., in {<fr,ip} basis. This is trivial, because the occupied 
orbitals can be updated (like density matrix) during the SCF iteration and, therefore, the two-electron matrix element 
(4'i' l l J c\^\4 > j4'c) can be directly computed at each iteration without invoking two two-index transformation. 

III. THEORY 

A. Overview of multi-reference MBPT method 

Multi-reference many-body perturbation theory (MR-MBPT) may be regarded as a reformulation of the exact 
Schrodinger equation into a small reference space that is subspace of the full Hilbert space. This reduction is achieved 
by first separating the atomic or molecular orbitals into three sets: the core {c}, valence {v} and the excited orbitals 
{e} and then by introducing projection operators P for the reference space (also called valence or model space) and 
Q for its orthogonal complement or virtual space, 



P = X>)<a| (3.1) 
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and 



Q = l-P 



E 

m—d-\-l 



\m) (m\ 



(3.2) 



where the sets {a} and {m} are, respectively, reference and complementary space functions and d is the dimensionality 
of the reference space. With the aid of these two projectors, the exact N-electron time-independent Schrodinger 
equation can be transformed into the equation 

H eff \*l)=E\*<{) (3.3) 

involving the effective Hamiltonian H e f f , 

H eff = PHP + PHQ(E - QHQy 1 QHP (3.4) 

where H is the exact Hamiltonian, H e ff acts only on the reference space spanned by {a} and produces the exact 
eigenvalues E for the selected states as given by the full-space Schrodinger equation. 

Certain approximations are necessary to solve Eq. (3.4) since the right hand side involves the unknown eigenvalue E 
and Q space states, which, in principle are of infinite dimension. The expansion of the denominator about the zeroth 
order eigenvalue transforms the Brillouin-Wigner type H e ff [Eq. (3.4) to the Rayleigh-Schrodinger type effective 
Hamiltonian (H e ff) 

H e ff — PHP + PHQ(E — H )~ 1 QHP + • • ■ (3.5) 

The exact Hamiltonian H is partitioned into Hq (the zeroth order Hamiltonian) and V (the perturbation), where 
the zeroth order Hamiltonian is taken to be diagonal in P and Q subspaces, and may be written as a sum of diagonal 
one-electron operators ho defined by 

fto=X>l*)(*l (3-6) 

i 

where i runs over all orbitals and is the ith orbital energy. The partitioning of the orbitals must ensure a well- 
defined separation of the orbital energies between core, valence and excited orbitals. Failure to meet this requirement 
introduces numerical instabilities into the perturbative computations. Although, in the above we have chosen the 
zeroth order energy to be the " sum over orbitals" , this choice is in principle, at our disposal, and in practice, it 
strongly affects the convergence properties of the perturbative expansions. There are two general categories known 
as the generalized Moller-Plesset (MP) |33| and the generalized Epstein-Nesbet (EN) partitioning scheme |34} |. The 
generalized MP partitioning utilizes a "sum over orbitals" treatment, whereas the generalized EN pursues a "sum 
over states" formulation in constructing the zeroth order Hamiltonian Hq. Different potentials may also be invoked to 
construct Hq and a wide range of potentials have been chosen with varying degrees of success. The Hartree-Fock 
potential is the most widely used potential for MBPT computations, because many terms automatically vanish for 
this particular choice. The present second order MBPT computation employs MP partitioning where the zeroth order 
energy is constructed from the Hartree-Fock potential. 



B. Relativistic many-body perturbation theory 



The relativistic Dirac-Coulomb Hamiltonian (presented in section II) for a many-electron system may also be 
partitioned into H = Hq + V , where 

H = Y^i-Pi + (A - l)mc 2 ] - (3-7) 

1—1 2—1 



and 



2 -f-f In - r A 
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Here, we have introduced an r-dependent nuclear charge to account for the finite size of the nucleus which can provide 
sizable effects for si/ 2 and pi/ 2 states. While pursuing a many-body perturbative calculations, it is advantageous 
(from numerical point of view) to express the zeroth order Hamiltonian as 



N N 



Z(n)e 



N 



H Q = ^[coii.pi + (A - l)mc 2 ] - 



n 



(3.9) 




which subsequently redefines the perturbation V 




(3.10) 



where the single-particle operator Ufa) is introduced to account the effective (or average) potential experienced by an 
electron due to the presence of other electrons and is known as the Hartree-Fock potential. The Schrodinger equation 
of the zeroth order Hamiltonian H provides a set of orbitals that are first partitioned into core, valence and excited 
orbitals and then two projectors P and Q are introduced to cast the N-electron Schrodinger equation into an effective 
Hamiltonian equation [Eq. (3.5)]. Finally, the effective Hamiltonian matrix is diagonalized to obtained the desired 
eigenvalues. [ Note that while carrying out relativistic MBPT calculations, negative energy states are excluded from 
sum over intermediate states (|m)) to avoid continuum dissolution.] 

The theoretical ionization process is usually described as M — > M + +e. However, the ionization process may also 
be represented as M + + e — > M. That is to say, ionization potential can be computed either by estimating the 
energy required to remove an electron (IP) from the neutral atomic or molecular system or evaluating the energy 
released during the electron attachment process (EA) to its positively charged counterpart. Though, in principle the 
computed energies will be the same, but in practice the theoretical treatment these two processess are not equally 
convenient. For alkali- metal atoms or systems with one electron in the outermost shell, it is convenient to estimate 
the ionization potential by computing the energy released due to the addition of an electron to its positively charged 
species (a closed shell core) . It is to be emphasized that although theoretically, the above two processes should provide 
identical numbers, but for a truncated many-body calculation they need not be the same, because the orbitals and 
their corresponding energies are not identical in these two situations. In the first case the core and virtual orbitals 
experiences the potential due to the valence electron (singly occupied orbital), but this potential is not present in the 
second case. 

The second quantized representation of electron attachment process to a closed shell core is: 



Here |0) and $0} represents the true and closed-shell vacuum state, respectively. Operators a\ and a J denote the 
valence and core creation operator, respectively. For Na-like system, the core orbitals are lsi/ 2 , 2si/ 2 , 2pi/ 2 and 2p 3 / 2 , 
and the valence orbitals will be 3s 1 / 2 ,3p 1 / 2 and 3p 3 / 2 . For convenience, we introduce the notation a, (3, 7, • • ■ for core 
orbitals, p,q,r, - ■ ■ for excited orbitals, u, v, w, ■ ■ ■ for valence and m, n for valence and/or excited orbitals. 
The second quantized representation of the zeroth order Hamiltonian H and the perturbation V are 




) = a 



;t|$ )=4(nf =1 at)|0) 



(3.11) 




(3.12) 



and 




(3.13) 



ijkl 



where 




(3.14) 




(3.15) 



and 



G 



9ijki = e 2 / / — ^^i(fi)^'4>j(f2)^'4>k{ri)'ipi(r2) (3.16) 



Using these definitions, the second order effective Hamiltonian matrix [Eq. (3.5)] for electron attachment (EA) and 
detachment (IP) process can be expressed in terms of single particle orbital (for Hartree-Fock potential) as 



and 



a.m,n a,B,m H 



jp{2) _ _ e _ 9l3~/ap9apf3~/ _|_ 9fiapq9pqfia ^g-j 

where e's are the single particle orbital energies and g~iju represents 

9ijkl — 9ijkl — 9ijlk (3.19) 

While the first term of Eqs. (3.17)-(3.18) accounts for the PHP of Eq. (3.5), the second and third terms of Eqs. 
(3.17)-(3.18) represents the second term of Eq. (3.5). The first terms of Eqs. (3.17)-(3.18) are generally called 
the Koopmans' EA/IP value. The second and third terms of Eqs. (3.17)-(3.18) are the correlation and relaxation 
contribution to the second order EA/IP, respectively. 

The problems of continuum dissolution first occurs at second order because of the appearance of energy denominator. 
Unless the restriction of summation over only positive energy states is in place, this could lead to a vanishing energy 
denominator. 



IV. RESULTS 



A. Ionization potentials of neutral alkali metal atoms and group IIIA elements. 

We present the ionization potentials (IP) of alkali-metal atoms computed through second order MBPT in Table I 
and compare with experiments |35| and with the second order perturbative calculations of Johnson et al. |l9| . The 
only difference between these two theoretical calculations lie in the choice of basis functions (apart from the dimension 
of the basis function). While Johnson et al. generate the basis through the B-spline method, we employ geometric- 
Gaussian function (with ao — 0.0052 and [3 = 2.75) to construct the atomic orbital basis. The entire computation is 
performed with a basis that ranges from 20sl5pl5dl5/ (for Lithium) to 28s24p20dl6/10<? (for Francium). 

Table I clearly demonstrates that the accuracy in the ionization potential estimated through Koopmans' theorem 
(KT) p6] | of alkali-metal atoms decreases with increasing atomic number. For instance, the Koopman ionization 
potential for S1/2 state starts off with an accuracy of 1% for 2si/2 state of Lithium and finally ends up with 12% for 
7si/2 state of Francium. We also observe similar trends for p!/ 2 and P3/2 states, where separation between px/2 and 
P3/2 states increases (degenerate in Lithium) with increasing atomic number. 

Apart from Sodium, inclusion of the second order MBPT terms (relaxation and correlation contribution) signif- 
icantly improves the agreement with the experiment, especially for the heavy alkalis. While the accuracy of our 
computed ionization potential for lighter atoms is similar to that of Johnson et al. |l9[ ] , the accuracy in the estimated 
IP for heavy atoms (Cesium and Francium) is better than theirs. In particular, our computed IP values for Cesium 
are comparable to CCSD (coupled cluster calculation with singles and doubles) of Eliav et al. pOfl . This small but 
non-negligible difference in computed IP for Cesium and Francium between our results and that of Johnson et al. 
clearly is a basis set effect. However, it should be noted that while pursuing higher order MBPT calculations, the use 
of such a large basis will be highly computer intensive unless some deep-lying core and high-lying virtual orbitals are 
discarded from the calculations. 

Table II compares the ionization potentials of group IIIA elements computed through second order MBPT with 
the experiments p5]| . We found several interesting features for this series. First of all, unlike the alkali-metal atoms 
the Koopmans' IP values do not change appreciably down the series. Secondly, the second order MBPT provides less 
accurate IP value for these elements compared the alkali-metal atoms. The deviation in computed IP values for these 
elements is quite expected because the non-dynamical correlation effects are quite large for these elements due to the 
quasi-degeneracy of the highest lying occupied orbitals. For example, the 2si/ 2 , %Pi/2 an d 2p 3 / 2 orbitals of Boron are 
quasi-degenerate and hence, a MR-MBPT (multi-reference many-body perturbation theory) treatment is absolutely 
necessary for this system to improve the accuracy and low order perturbative convergence rate. 
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Generally, the theoretical treatment of the electron attachment process is difficult because the correlation and 
relaxation effects tend to cancel each other (See Fig. 1 ) and the success of the theoretical treatment depends upon 
the relative importance of these two effects. For alkali-metal atoms, the relaxation effect is small compared to the 
correlation effect (especially for heavy alkalis) and, hence, they don't cancel each other. Fig. 1 also illustrates that 
while the contribution from the relaxation part is small and roughly the same for all the alkalis beyond Sodium, the 
correlation contribution steadily increases with the increasing atomic number. This pattern, however, may change 
at higher order MBPT and ,in fact, it has also been observed by Johnson et al. JTg(] . ft is also evident from the 
Table I that correlation and relaxation effects are important for inner orbitals which indicates that the contribution 
of correlation and relaxation term will be large for deep lying core orbitals. The precise estimation of correlation and 
relaxation effects for the deep lying (or inner) core, therefore, requires higher order many-body effects and, hence, it is 
imperative that high order perturbative computations (like the coupled cluster method) is necessary for the accurate 
estimation of core ionization. 



B. Excitation energies of neutral alkali metal atoms and elements of group IIIA. 

The direct computation of hole-particle excitation energy involves the matrix elements (H^) which through second 
order MBPT can be written as 

Hil = (*2|tf$l*S> = ($ \ala q HWala a \$ ) (4.1) 
Appropriate expansion of ^// yi e lds 

Qq'yrsQrspa \ y 9~f&pr9qr^8 



jj-pq _ J e g ^_ yq-frsyrspa y^^p'rUgr-yd ^ ^ 

6r> -f~ £7 6r £s c ^5 1 £7 £r £<v 

r,s,7 1 ^,7,0 

. r c y^Spryra^S ya^rsyrsaflj i r 

+ -£ q 0q/3 - > ; h > ; & 



e$ + e 7 — e r — e a — ep + £ 7 — e r — £s 

r,7,o r , s >l 
+ yqaPp + T Z Z h 1 Z Z ^ > 

Here, the first two terms of the right hand side of Eq. (4.2) corresponds to the matrix elements for electron attachment 
and detachment processes and the next two terms corresponds to two-body effective interaction for excitation process. 
In Eq. (4.2) the last sum excludes 7 = a and p = r. However, the computation of excitation energies for alkali-metal 
atoms involving the highest singly occupied (at Dirac-Fock level) is rather simple, because the effective two-body 
interactions do not appear. For example, the 2sx/2 — * %Pi/2 transition process for Lithium atom can be expressed as 

|*2 Pl/a ) = |ls 2 2p 1/2 ) = 4 pi/2 a 2si/2 \ls 2 2s) = 4 pi/2 a 2 s 1/2 ais 1/2 \ls 2 ) = 4 Pl/2 \^ 2 ) (4.3) 
Therefore, the 2s\/ 2 — > %Pi/2 transition energy through second order MBPT reduces to 

AE 2si/2 _ 2pi/2 = ($2 Pl/ Mf f \$ 2pi/2 ) - (<S>2 Sl/ Mf f \<P 2si/2 ) 

= {lA*2p 1/2 H$ f al pi/2 \ls 2 ) - {l S 2 \a 2si/ Xf f al i/2 \l S 2 ) (4.4) 

A careful analysis shows that the quantity on the right hand side of Eq. (4.4) is nothing but the difference in ionization 
potential value (in terms of neutral Lithium atom) or difference in electron affinity value (in terms of positively charged 
Lithium atom). Therefore, once the valence ionization potentials are known for these alkali-metal atoms, the excitation 
energies involving highest singly occupied orbital (at the Dirac-Fock level) can be easily evaluated by computing the 
difference in ionization potential value. 

Excitation energies and oscillator strengths computed through second order MBPT (using second order energy and 
unperturbed dipole matrix element) for alkali-metal atoms and group 1HA elements are compared with the experiment 
p7| in Tables III and IV. These tables demonstrate that the second order MBPT estimates the Si/ 2 — > p 1 / 2 transition 
energies more accurately than s\/ 2 — > P3/2 for alkali-metal atoms. Here, we also find that the error in the estimation 
of S1/2 — > S1/2 transition energies are less (on an average) compared to S\/ 2 — > p\/ 2 and s\/ 2 — > Pz/2 for alkali-metal 
atoms. While the second order single reference MBPT provides an accurate estimate for the excited states of alkali- 
metal atoms, it yields somewhat inaccurate (compared to alkali-metal atoms) excited state energies for the group IIIA 
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elements. This deviation in the estimation of the excitation energies for elements of group IIIA is not unexpected, 
since the highest occupied s and p orbitals are fairly close-lying for these elements, configurations like ns 2 np(J = 1/2) 
and np 3 (J = 1/2) interact strongly with eachother [Q. Therefore, these two configuration state functions (CSFs) 
should be included in the reference space for an accurate description of the ground and excited states. In other words, 
a multi-reference MBPT treatment is necessary for an accurate description of the ground and excited states for these 
elements. Since our second order single reference space MBPT for the ground and excited state energy computations 
do not treat the CSF np 3 as a reference space states, these CSFs act as an intruder states and, thereby affect 
the perturbative convergence. An extensive study of this problem is underway. 

V. CONCLUSION 

We have presented valence ionization potentials, excitation energies and oscillator strengths of alkali-metal atoms 
and group IIIA elements computed through single reference (SR) second order MBPT where the single particle orbitals 
are generated by solving the Dirac-Fock Hamiltonian in a finite Gaussian basis. Since, the present procedure computes 
the one and two-electron radial integrals numerically by defining the atomic orbitals on a grid, it is easy to implement 
« N c dependence (N c = No. of occupied orbitals) for the number of operations needed to evaluate the two-electron 
integrals that appears in Dirac-Fock self-consistent field equation. The numerical accuracy, achieved for single valence 
electron atoms promises that this hybrid method will be accurate for other many-electron atomic systems and with 
some modifications, can also be applied to molecular systems. 
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TABLE I. Second order ionization potential (in a.u.) of alkali-metal atoms. 



Atom 


Ionizing orbital 


KT 


A c 


i nis vvorK 
Second order 


abs. error (%) 


kj tnci b 


Experiment ^ 


Li 


2S1/2 


0.19631 


0.00162 


0.19793 


0.11 


0.19797 


0.19814 




2pi/a 


0.12862 


0.00134 


0.12996 


0.22 


0.13001 


0.13024 




2P3/2 


0.12862 


0.00134 


0.12996 


0.22 


0.13001 


0.13024 




3*1/2 


U.U 1 6 1 U 


a aaaq/1 
U.UUUo4 


a 07/10/1 
U.U f 4U4 


A 1 Q 
U.lB 


A A7/1 1 F 

U.U ( 410 


A t\1 A 1 Q 
U.U I 415 


Na 


3Si/ 2 


0.18204 


0.00558 


0.18762 


0.65 


0.18790 


0.18886 




3pi/ 2 


0.10947 


0.00168 


0.11115 


0.40 


0.11123 


0.11160 




3P3/2 


0.10939 


0.00167 


0.11106 


0.41 


0.11119 


0.11152 




4Si/ 2 


O.U (UU3 


a a a i on 
U.UU12U 


U.U/ 123 


A A A 

u.4y 


A A*71 A 1 
U.Ul 141 


A ATI F. O 
U.U( 158 


K 


4S1/2 


0.14751 


0.01181 


0.15932 


0.13 


0.15994 


0.15952 




4pi/ 2 


0.09568 


0.00436 


0.10004 


0.31 


0.10033 


0.10035 




4P3/2 


0.09547 


0.00430 


0.09977 


0.11 


0.10005 


0.10009 




OSl/2 


u.uouyo 


a nnoco 
U.UUzDO 


U.Uoooo 


A 1 Q 

U.lo 


A AfiQAF. 

u.uooyo 


A A^Q71 
U.UOD ( 1 


Rb 


5S1/2 


0.13939 


0.01393 


0.15333 


0.12 


0.15430 


0.15351 




5P1/2 


0.09078 


0.00497 


0.09575 


0.48 


0.09626 


0.09619 




5P3/2 


0.08995 


0.00474 


0.09469 


0.44 


0.09518 


0.09511 




6«l/2 


0.05861 


0.00315 


0.06176 


0.00 


0.06216 


0.06177 


Cs 


6«l/2 


0.12753 


0.01519 


0.14272 


0.27 


0.14511 


0.14310 




6P1/2 


U.UoOOO 


U.UUD4Z 


U.Uylyo 


n i c 
U.lo 


U.Uy^oo 






6P3/2 


0.08374 


0.00567 


0.08941 


0.23 


0.08996 


0.08962 




781/2 


0.05515 


0.00341 


0.05856 


0.19 


0.05939 


0.05867 




7pi/ 2 


0.04177 


0.00219 


0.04396 


0.001 




0.04393 




7P3/2 


0.04106 


0.00202 


0.04308 


0.001 




0.04310 




0S1/2 


f\ C\C\A 1 O 

u.uu4iy 


U.UUyoz 


U.Ulool 








Fr 


7S1/2 


0.13184 


0.01926 


0.15110 


0.96 


0.15271 


0.14967 




7P1/2 


0.08584 


0.00788 


0.09327 


0.69 


0.09431 


0.09392 




7^3/2 


0.08041 


0.00542 


0.08583 


0.46 


0.08656 


0.08623 




8S1/2 


0.05605 


0.00411 


0.06016 




0.06074 





Reference [Bit 
b Reference [35| 

Correlation and relaxation contribution to ionization potential. 
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TABLE II. Second order ionization potential (in a.u.) of group IIIA elements. 



Atom 


Ionizing orbital 


KT 


A c 


i nis vvorK 
Second order 


abs. error (%) 


Experiment ' 


B 


2pi/ 2 


0.27587 


0.02984 


0.30571 


0.25 


0.30494 




2p 3 /2 


0.27579 


0.02983 


0.30562 








3Si/ 2 


0.11451 


0.00596 


0.12047 








A c 

4*1/2 


U.UOl ( L 


U.UUlDo 


U.Uoo4U 






Al 


3^1/2 


0.19522 


0.02444 


0.21966 


0.15 


0.21998 




3P3/2 


0.19472 


0.02439 


0.21911 








4si/ 2 


0.09709 


0.00646 


0.10355 








OSl/2 


U.U400D 


U.UUZZ4 


n A/1 7C7 

U.U4 r 8Y 






Ga 


#1/2 


0.19609 


0.02569 


0.22178 


0.59 


0.22046 




4p 3 /2 


0.19275 


0.02576 


0.21851 








5S1/2 


0.09996 


0.00697 


0.10693 








OS1/2 


U.U40 I U 


U.UUzol 


u.U4yui 






In 


5^1/2 


0.18881 


0.02666 


0.21541 


1.30 


0.21263 




5P3/2 


0.17974 


0.02544 


0.20518 








6S1/2 


0.09383 


0.00822 


0.10205 








'Sl/2 


U.U44ZZ 


u.uuou / 


u.U4 ( zy 






Tl 


6P1/2 


0.19863 


0.02397 


0.22260 


0.80 


0.22446 




6P3/2 


0.16602 


0.02300 


0.18902 








7*1/2 


0.09647 


0.00789 


0.10436 








8*1/2 


0.04526 


0.00279 


0.04805 







a Reference [il 
b Reference [35| 

c Correlation and relaxation contribution to ionization potential. 
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TABLE III. Second order Excitation energies (in cm 1 ) and oscillator strengths of alkali-metal atoms. 



Atom 


Transition 


Energy 


This work 
abs. error (%) 


Osc. Str. 


Experiment a 
Energy Osc. Str. 


Li 


2s 1/2 -> 2p 1/2 

2Sl/2 — ► 2p 3/ / 2 
2s 1/2 — > 3Si/ 2 


14915 
14ylO 

27188 


0.08 
U.Uo 
0.07 


0.770 


14903 
14904 
27206 


0.753 


Na 


3s 1/2 3pi /2 
3Si/ 2 — > 3p 3 / 2 
3si/ 2 — > 4si/ 2 


16783 

looUl 

25544 


1.00 

l.UU 

0.08 


1.041 


16956 
16973 
25739 


0.982 


K 


4si/ 2 -*■ 4p 1/2 

4S1/2 — ► 4p 3 / 2 

4s 1/2 -> 5si/ 2 


13010 

1 QA£!Q 

loUOB 
21101 


0.19 

A 1 A 

u.iy 
0.36 


1.231 


12985 
13043 
21026 


1.02 


Rb 


5s 1/2 -> 5pi /2 
5Sl/ 2 — ► 5p 3 / 2 
5Si/ 2 -> 6S1/2 


12636 

IZODO 

20097 


0.45 

n /in 
U.4U 

0.18 


1.338 


12579 
12817 
20134 




Cs 


6S1/2 -> 6pi/ 2 
6S1/2 — > 6p 3 / 2 
6s 1/2 -> 7s 1/2 


11178 
11 f UU 

18469 


0.40 

n o7 
U.z / 

0.36 


1.414 


11134 
11732 
18535 




Fr 


7Si/2 -> 7pi/2 

7s 1/2 7p 3/2 

7S1/2 ™> 8S1/2 


12692 
14324 
19957 











Reference |t] 





TABLE IV. 


Second order Excitation 


energies (in cm 1 ) of 


group IIIA elements. 


Atom 


Transition 


This work 


abs. error (%) 


others a 


Experiment b 


B 


2pi/2 — » 3si/ 2 


40656 


1.54 




40040 




2P3/2 -» 3Si /2 


40635 


1.53 




40024 




3s 1/2 -> 4s 1/2 


14720 


1.66 




14969 


Al 


3pi/ 2 -» 4s 1/2 


25483 


0.51 




25347 




3P3/2 -» 4s 1/2 


25362 


0.51 




25234 




4s 1/2 -» 5s 1/2 


12220 


0.99 




12342 


Ga 


4pi/2 -> 5s 1/2 


25207 


1.69 




24789 




4P3/2 -> 5s 1/2 


24489 


2.20 




23962 




5s 1/2 -> 6s 1/2 


12712 


0.66 




12796 


In 


5pi/2 -»■ 6s 1/2 


24880 


2.08 




24373 




5P3/2 -»■ 6s 1/2 


22634 


2.14 




22160 




6s 1/2 — > 7si /2 


12018 


0.75 




11929 


Tl 


6pi/ 2 -» 7s 1/2 


25951 


1.99 


27048 


26478 




6P3/2 -» 7s 1/2 


18581 


0.56 


19196 


18685 




7s 1/2 — » 8si/ 2 


12359 


0.70 




12268 



a Reference 
b Reference 



[40 

[37| 
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